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ABSTRACT 

This is a preliminary study to examine the prospect of detecting TeV photons from 7-ray 
bursts (GRB) using km-size neutrino telescopes, specifically for the ANTARES neutrino tele- 
scope. Although optimized to detect upgoing neutrino-induced muons, neutrino telescopes 
nevertheless have a potential to detect high-energy photons by detecting downgoing muons 
from the electromagnetic cascade induced by the interaction of TeV photons with the Earth's 
atmosphere. The photon energy spectrum of a GRB is modeled by a simple power law and 
is normalized by simple energy considerations. Taking into account the absorption of TeV 
photons by cosmic infrared backgrounds, an optical depth table calculated from a model by 
[Finke, Razzaque & Dermer|(2010) is used and the arriving number of photons on top of the 
Earth atmosphere is determined. Muon production in the atmosphere is determined by con- 
sidering two main channels of muon production: Pion photoproduction and direct muon pair 
production. The muon energy loss during their traverse from the surface to the bottom of the 
sea is determined using the standard muon energy loss formula. Assuming different detector 
sizes, the number of detectable muons from single GRB events located at different redshifts 
and zenith distances is determined. The background is calculated assuming it consists primar- 
ily of cosmic ray-induced downgoing muons. The detection significance is calculated and it 
can be concluded that to obtain at least 3cr detection significance, a typical GRB has to be 
located at redshift z < 0.07 if the detector's muon effective area is A^^ ^ 10^^ km^, or 
redshift z < 0.15, if the muon effective area is A^'jj ^ 1 km^. 

Key words: gamma-ray burst: general — elementary particles — nuclear reactions, nucle- 
osynthesis, abundances — methods: analytical 



1 INTRODUCTION 

A downgoing TeV 7-ray passing through the Earth's atmosphere 
will initiate an electromagnetic shower that contains a small num- 
ber of muons. A kilometer- size neutrino telescope can detect these 
muons by detecting the Cherenkov photons produced from the in- 
teraction of these downgoing muons with the medium surrounding 
the telescope. The signals can then be reconstructed in order to ob- 
tain the photon's energy and its direction of arrival. This method of 
detection can, in principle, be applied to various TeV 7-ray sources 
but for reasons that will be discussed in the following, this paper 
concentrates on detecting TeV photons from 7-ray bursts (GRB). 
This paper is part of the theme of GRB detection by neutrino tele- 
scopes. Neutrino-induced muon event rates has been calculated 



in detail (e.g.lHalzen & Hoo per|1999 Alvarez-Muniz, Halzen &| 
|Hooper|2000l[Gupta|2002(|Dermfi'r & Atoyan|2003^ , and this paper 
is complementary to them by providing the necessary calculations 
of photon-induced muon event rates. 



Since the first detection of 7-ray burst (GRB) of cosmic ori- 
gin ( [Klebesadel, Strong & 01sonl|1973) in 1969, many advances 
has been made contributing to the study of GRBs (see Fishman 



1995';'Fish man & Meegan| 1995 [[van Paradijs, Kouveliotou & Wi 
fjers 2000 for excellent reviews of the early days of GRB astron- 
omy). Many challenges have been answered and problems solved, 
but many questions still remain and one of the outstanding prob- 
lems in TeV 7-ray astronomy is the detection of TeV photons emit- 
ted from GRBs. Many authors have constructed models that pre- 



dict the existence of TeV emission (e.g. Meszaros, Rees & Pap' 
athanassiou 1994, Dermer, Chiang & Mitman 2000, Pe'er & WaX' 
|man|2005(,Asano & Inoue,2007) , and there are no indications of a 
cut-off energy ever observed in any GRB spectrum. It is therefore 
possible that TeV photons could be produced and if detected they 
can provide constraints and clues of the GRB production mecha- 
nism. Not only their detection can give clues on the intrinsic prop- 
erties of GRBs, but they can also provide constraints on the ex- 
tragalactic infrared background that attenuate TeV photons during 
their travel from the source to the observers. 
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Observational evidence for higher-energy photons has accu- 
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mulated in the past years. BATSl|^ observed several thousands 
GRBs with photon energies up to 300 keV ( [Paciesas et al.|1999| l, 
followed later on by EGRElj^(an instrument on board the Compton 
Gamma Ray ObservatorTj observations of GRBs with energies up 
to 200 MeV ( [Dingus|199"5i|Schneid et al.|1992^ and even a peculiar 
GRB that la sts for 90 minutes and energies up to 18 GeV jHurley| 
et al.|l994 i. Most recently, FermF] observed GRBs with energies 
up to ~ 10 GeV fAbd o et al.|2009t . 

Moving higher up the energy scale, attempts have been made 
to detect TeV components of GRBs. Using coordinates distributed 
by the BATSE Coordinates Distribution Network (BACODINE) 
and later on by the GRB Coordinates Network (GCN), the Whip- 
p/fQcoUaboration has observed 9 BATSE GRBs and 7 other GRBs 
announced by GCN within minutes to hours after the burst time 
given by the alert ( [Connaughton et aL]|1997[ [Horan et al.||2007| l. 
No evidence of TeV emission were found but upper limits were re- 
ported. The MAGIC^Telescope, using the same observation prin- 
ciple as Whipple, observed 9 GRBs announced by GCN and found 
no evidence of TeV emission as well ^Albert et al.|2007l l. 

So far the only indication of TeV emission were detections 
by the HEGR^ AIROBICC fl detector who claim ed > 16 TeV 
emission from GRB 920925clPadilla et al.|l9"98t, the Milagritcj^ 
collaboration ^Atkins etal . 2000 2003 , 20051 who report d etection 
at 0.1 TeV, and the GRANEj^ array (Poirier et al.|2003| l at 0.01 
TeV. These reports, however, have marginal detection significance 
(~ 3(t) but nevertheless provide a tantalizing glimpse of the possi- 
ble extension of the GRB energy spectrum. 

Undersea or under-ice neutrino telescopes are far less sensitive 
to TeV photons than Imaging Atmospheric Cherenkov Telescopes 
(lACTs). But why it is argued that km-size Cherenkov detectors 
such as ANTARE:^ IceCub{(3 or the future KMSNelj^have 
the potential to observe TeV photons from GRBs is not only be- 
cause of its large collecting area but also because of its wide field 
of view and high duty cycle. Its field of view can reach up to ~ 3 sr 
above the horizon and operate at ~ 95% duty factor. This means 
that the detector is taking data most of the time and will almost 
always be able to record events from any nearby GRB above the 
horizon which emits photons numerous enough to induce muons in 
sufficient number to be detectable. On the other hand, lACTs can 
only operate on clear, moonless nights, and their slow slewing rate 
is inadequate to observe GRBs which are transient by nature. More- 
over, if we know the time and direction of the burst, we then have 
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the knowledge of when and where to observe and this can greatly 
reduce the amount of background of downgoing muons induced by 
cosmic rays. This two-piece of information can be provided by the 
GRB Coordinates Network (GCnQ a system that distribute alert 
notices to its subscribers whenever any spacecraft that is part of this 
network detect a potential GRB l |Barthelmy et al.|2000t . 

At present fiv e satellites ar e part of this network: HETE ^ 
(Ricker et al.'2003'), INTEGRAlPj ( Winkler et al. 2003), Swi^ 
l ,Gehrels et al..2004j , Fermi ( .Moiseev 2008 ), and AGIL^lCocco 



|et al.|[2002l l. By performing time and position correlation of de- 
tected events with those provided by GCN we can significantly re- 
duce the background and increase the possibility of detection at 
higher significance. 

Despite these potentials, detecting the TeV component of 7- 
ray bursts is not without pitfalls. One of the main problems that 
comes to mind is the attenuation of TeV photons by ambient IR 
photons in the universe. Along their path from the source to the 
Earth, TeV photons collide with ambient IR photons and annihi- 
lating themselves, creating pairs of electron-positron. The cross 
section for such process is well-known but measuring the accu- 
rate spectral density of cosmic IR photons at all redshifts is still the 
main problem. I will discuss this problem in more details in Section 
|2.2| by confronting current attenuation models with observational 
data and choose the best model. This attenuation will consequently 
limit our observations only to the nearest GRBs. 

Another crucial problem is to calculate the number of muons 
produced from a 7-shower with a certain spectral index. Two muon 
production mechanisms are identified: photoproduction and di- 
rect muon-pair production. Both mechanisms are low-cross-section 
process and are efficient at different energy regimes. Both mecha- 
nisms will be discussed in more detail in Section[3]and I will pro- 
vide the necessary formula to determine the number of muon pro- 
duced from 7-showers. In calculating the observed muon flux at de- 
tector level, I also take into account the muon energy loss caused by 
their passage through seawater (Section|4|. Using all this, the num- 
ber of detectable muons for single GRB events at different redshifts 
are calculated. This is outlined in Section|5] 

Previous works had addressed these problems. In a similar 
vein with this paper, |Alvarez-Mufiiz & Halzen| ( [1999 ') have cal- 
culated the muon flux from GRBs at various redshifts and phys- 
ical properties, but for the muon production only the photopion 
channel — whic h falls quickly with increasing energies — i s consid- 



ered. Recently, Halzen, Kappes & 6 Murchadha ^2009| rectified 



this problem and discussed the various channels from which muons 
can be produced, including muon pair production. However they 
only consider nearby, local, TeV sources and neglect altogether 
extragalactic sources. Furthermore their results indicate that local 
TeV sources are too faint, even an IceCube-size neutrino telescope 
require ~ 10 years of integration time before a detection with Scr 
significance is obtained. 

This paper is the first part of a series of papers that try to ad- 
dress these problems and examine the possibilities of detecting TeV 
photons from GRBs. This paper (Paper I) will outline the processes 
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that take place during the emission of TeV photons up to the de- 
tection of muons at detector level and will also formulate the basic 
working equations. Paper II will present a calculation of the muon 
event rate from stacked GRB sources, made using the latest distri- 
bution functions of GRB world model, and Paper III will present a 
data analysis from the ANTARES neutrino telescope. 



2 THE PHOTON SPECTRUM OF A GRB 

2.1 Normalization of tlie photon spectrum 

The photon spectrum of a GRB occuring at redshift z is assumed 
to be constant during the whole duration of the burst. The burst du- 
ration in the observer's frame is At = (1 + 2) At, (throughout this 
paper, asterisks will be used to indicate terms in the source's frame 
while terms without asterisk are terms in the observer's frame). 
The photon spectrum N{e) of a GRB is approximated by a broken 
but smoothly connected power law, known as the Band spectrum, 
which is a model based on BATSE observations of 54 GRB ijBandj 
|etal.|I993| >: 



V^bk 



H{e - Ebk) exp(a -h)[ — 

\ Ebk 



-(b+1) 



-(a + 1) 

TeV^^ cm^^ s" 
(1) 



where (a, 6) are respectively the spectral indices of the power law 
in the low- and high-energy regime demarcated by the break energy 
ebk, and is the normalization constant in unit of photons TeV~^ 
cm~^ s~^. The function H{x) is the Heaviside step function de- 
fined as H{x) = 1 for a; and H{x) — otherwise. 

The break energy is related to the directly measurable peak 
energy tpk, , which is the energy in which the uf^ = e^N{(.) peaks, 
through 



Ebk = 



-Epk- 



(2) 



BATSE observations extend only to several hundreds keV and 
in some cases to several MeV, but subsequent observations by later 
satellites confirmed that the power law also extends to at least sev- 
eral GeV (e.g. |Hurley et al.| 1 994[ [Gonzalez et al.|2003||Abdo et al.| 
|2009[ >. Based on this observational evidence, in this work we as- 
sume that this power law function also extends to the TeV regime. 

The normalization constant is calculated by relating the 
energy spectrum in Equation[T|to its instrinsic isotropic-equivalent 
bolometric luminosity -Lboi, : 



(3) 



Lt°u = A-Krl{z){l + z) dt deN{e)e. 

Jo Jo 

in which rc{z) is its comoving distance at redshift z: 

r.iz)^£dz'il^^')^, 

where dl / dz is the cosmological line element defined as 

dl c 1 



dz 



Ho (1 + z)y^nZ+aJYTW' 



(4) 



(5) 



in which c is the speed of light. Ho = 72kms^^ Mpc^^ is 
the Hubble constant at the present epoch, and (SlA,f7m) ~ 
(0.742,0.258) is the present dark energy and matter density in 



the universe in units of the critical energy density, 3H^ /SttG. It 
is assumed that the GRB emission spectrum is constant during the 
whole burst duration. It is also important to note that I/b°i* is an 
isotropic-equivalent luminosity which assume that the 7-ray emis- 
sion is isotropic and is not beamed. The true, beamed, bolometric 
luminosity L^^oil is related to Lboi, by 



itruc /-| 
bol* — U " 



„ n \ T true 

cos6^j)Lboi*, 



(6) 



where 6j is the opening angle of the jet. The average value of 



the opening ang le is (Oj) ~ 6° iGhirlanda erar]|2007[ >, making 



UoT. ~ o.oossLl:,'^"" 



The integration in Equation|3]is solved by locking the spectral 
index a to the typical value of a = ( [Preece et al.|2000]|Natarajan| 
|et al.|20()5} and letting the other values as free parameters. Solving 
the integration this way, we can obtain the photon flux f.y : 



(7) 



47rr2(z)Af,e2j^,Abor 

in which ebk* ~ ebk(l + z) is the break energy in the source's 
frame and Aboi is a bolometric correction to the flux, which is the 
result of the integration in energy. To avoid a divergent flux in the 
integration, we do not integrate it to infinite energy but instead cut 
the spectrum off at maximum energy emax* = 300 TeV. At the 
moment the upper cutoff of the photon spectrum is not known, and 
in fact the taking of 300 TeV as the limit of the integration is quite 
arbitrary. Taking this in mind, the value of Aboi is then 



Aboi 



-iexp 



cxp( — b) 
l-b 



(-fe) + |+exp(-fe)ln(^J^) 



for 6 / 1 

forfo = 1. 
(8) 



Thus given (-£/boi» 1 ^1 1 ^ bk» ) as parameters, we can construct 
the photon spectrum of any GRB. 

2.2 Photon absorption by ambient IR-photons 

Along the path from the source to the Earth, 7-ray photons inter- 
act with extragalactic background light (EBL) through the 77 — > 
e+e^ process, annihilating themselves and creating electron- 
positron pairs. A 7-ray photon with energy can produce a pair of 
electron-positron if it impacts a background photon with threshold 
energy 

2e^ 

eth = (9) 



where \ii — cos 9i is the angle of impact between the two photons. 
For head-on collisions, the wavelength of EBL photons which will 
interact with passing TeV photons is then 

Aebl ^ Ae— = 2.4e^ /^m, (10) 

in which Ae — h/{mec) is the Compton wavelength for an electron 
and the input e-y is in TeV. We can see that TeV photons will interact 
strongly with IR photons in the EBL. 

The optical depth r-y-, (e^ , 2) as a function of observed pho- 
ton energy and redshift z can be calculated if we also know the 
differential number density of background photons n(ebg, z) at en- 
ergy Ebg and redshift z: 



dz 



dl 
dz 



dfii(l — fii) 



(11) 



debgn(ebg, z)a[€^{l + z), tbg, Mi]. 
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in which emin — eth(l + 2)" , ^ is the cosmological hne ele- 
ment defined in Equatio n[5| and <j{ei, £2, fJ-i) is the cross section of 
the 77 pair production ( jBreit & Wheeler|1934[ |Gould & Sdireder] 
[T967|: 



(3-nin 



1-/? 



+ 2/3(/3^ - 2) 



(12) 



in which ctt is the Thomson cross-section and /3 is the electron- 
positron velocity in the center-of-mass frame: 



/3 



2^ 

eie2(l — 



^J■^ 



(13) 



Directly observing EBL photons to obtain their photon dis- 
tribution is difficult because of contamination issue from the in- 
strument as well as from the zodiacal light. Source discrimination 
is also another issue: The Cosmic Infrared Background (CIB) — 
which is extragalactic in nature — must be discriminated from fore- 
ground objects such as discrete sources like stars and compact ob- 
jects within the Galaxy, as well as diffuse sources such as light 
scattered and emitted by interplanetary dust and emission by inter- 
stellar dust (see |Hauser & Dwek| ( |2001^ for a review on this matter). 

There are many approaches in calculating the EBL photon 
density for all redshifts. One basic approach of doing it is by us- 
ing "backward models," in which we start from the existing galaxy 
count data and then model the luminosity evolution of these galax- 
ies backward in time (e.g. [Stecker, Malkan & Scully||2006) . An- 
other approach is the "forward evolution," performed by assuming 
a set of cosmological theory and semi-analytic merger-tree models 
of galaxy formation to determine the star formation history of the 
universe (e.g. |Primack, Bullock & Somerville|2005[[Gilmore et al.| 
|2009| l. Yet another approach is to focus on the properties and evo- 
lution of starlight, the primary source of CIB emission. This model 
integrates stellar formation rates and properties over time to obtain 
the amount of light emitted (e.g. jKneiske et al.|2004[|Finke, Raz^ 
|zaque & Dermer|2010) . 

In this paper I consider three different attenuation models: The 
"best-fit" model of Knei ske et al.| 12004 1, the fiducial model of 
Gilmore et al. (2009), and the recent "Model C" by |Finke, Raz-| 



zaque & Dermerj ||2010|). These mod els, along with the Baseline 



Model of Stecker, Malkan & Scully 1 2006 1, are compared in the 
plot of the Fazio-Stecker relation 1 Fazio & Stecker| 19701 in Fig- 



ure [T| The Fazio-Stecker relation is the (£7,2) value that gives 
r-y-y = 1. This is interpreted to be the redshift at which the flux 
of photons of a given energy is attenuated by a factor e and is 
called the 7-ray horizon. In this plot, for all models except those 
of [Stecker, Malkan & Scu lly (2006), for redshift < 5 the universe 
is optically thin to photons with energy < 20 GeV. At very low 
redshifts however, the models are relatively consistent with each 
other, but the differences start to become apparent at 2: > 1. The 
model calculated by Stecker, Malkan & Scully ( 2006) , which pre- 
dicts higher attenuation at higher redshifts, has in recent times con- 
tradicted MAGIC {Albert et al.|2008t and Fermi ( |Abdo et al.|20T0) 
observations and thus can be ruled out with high confidence (fur- 
thermore. Figure 2 in'Abdo et al. (2010) indicate that models by 
Finke, Razzaque & Dermer (2010); Gilmore et aL| ( |2009^ ; |Frances^ 
chini, Rodighiero & Vaccari (2 008] ) are the favorable ones) and will 



not be used in further calculations. 

Thus, knowing the attenuation function, we can then estimate 
the total number of photons emitted from a GRB at redshift z per 
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Figure 1. A plot of the Fazio-Stecker Relationship jFazio & Stecker|1970) 
for several attenuation models, as a function of redshift. Also shown are 
the redshifts and highest energy photons emax of various objects observed 
by Atmospheric Cherenkov Telescopes and Fermi-hAI jpinke & Razzaque] 
[2009, Abdo et al.|2010> . 
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Figure 2. An illustration of the effect of attenuation to a photon spectrum. 
Attenuation is calculated using the model by Fin ke, Razzaque & Dermer] 
(20101. The shape of the photon spectrum of a source located at redshifts 
indicated beside each curve is shown. Energies are in the observer frame of 
reference. The further a source is located, more attenuation is suffered by 
the highest energy photon. The curve is normalized to an arbitrary unit. 



unit energy arriving at the top of the Earth's atmosphere per unit 
area per unit time to be 



70(^7) = 1(^1,* = 0) = f-, 



-(6+1) 



(14) 



where is as derived in Equation [s] and 7(e, t) is the notation for 
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the photon flux at slant depth t in the atmosphere, as introduced in 
[Rossi & GreisenH1941^ . Slant depth t = Q means the top of the 
atmosphere. In this equation only the high-energy part of Equation 
[T]is used, because this is precisely the concern of this study and 
henceforth this equation will be the working equation. 

To give an illustration of the effect of attenuation to a pho- 
ton spectrum, the shape of the photon spectrum curve of several 
sources emitting at different redshifts is showed in Figure |2] As a 
comparison an unattenuated photon spectrum is also shown. The 
curves are normalized to an arbitrary unit. From the shape of the 
curves, the more distant the source is located, the more the photon 
spectrum curve is distorted due to attenuation effects. This imposes 
a limit on the number of TeV photons that we can observe. 



2.3 The cascade equation: Approximation A 

High energy photons interact with atoms in the atmosphere and ini- 
tiate electromagnetic showers of particles that will cascade on their 
way through the atmosphere. Through materialization or Compton 
collision, pairs of electron-positron will be produced, which in turn 
emit additional photons by way of bremsstrahlung. At each step the 
number of particles increases but their average energy decreases 
( [Rossi & Greisen[1941| l. Nevertheless these secondary photons can 
also produce muons that can be detected by the detector array, and 
thus it is important to calculate the total number of photons pro- 
duced in such a photon shower. 

This problem of counting particles produced in electromag- 
netic showers can be solved if we consider only radiation phenom- 
ena and electron-pair production, which can be described by the 
asymptotic formula for complete screening. This solution is called 
Approximation A ( [Rossi & Greisen|1941^ and allows us to calcu- 
late the photon flux at some depth t in the atmosphere, given the 
initial photon energy spectnmi. If the initial spectrum is in the form 
of a power law suc h as 'yje) oc e~^''"^^\ then the resulting spec 



trum at depth t is I Rossi & Greisen||l941 
[Murchadha|2009t 



Halzen, Kappes & O 



7(e,t) =7(e,t = 0) 



(ao + Ai)((To + M) 



oxp(Ait) 



A2 — Ai 
exp(A2t) 



OQ + A2 



(15) 



In this Equation as well as the in the following calculations, t 
is the slant depth in units of radiation length (in the atmosphere, 1 
radiation length equals 36.62 g cm~'^), cro = 7/9 is the probability 
per radiation length that an electron pair production will take place 
(in a case of complete screening), and A1.2 are the scale lengths 
factor of the shower growth and dissipation in the atmosphere. The 
formula to calculate A1.2 as a function of spectral index 6, as well 
as its tabulation, is given in [Rossi & Greisen| ( [1941 ). For 6 < 1, 
Ai is positive while for & > 1, Ai is negative. This would mean 
that in the former case the shower would grow as it penetrates the 
atmosphere while in the latter it will dissipate. Thus for a general 
case of an arbitrary value of b, the photon flux can be decomposed 
into its spectrum at the top of the atmosphere and its scale factor at 
depth t, i.e 



7(e,t) = 7o(<;)72(i)- 



(16) 



Particularly important is the case for fe = 1 since Ai = 
and A2 < 0, and this would make the second exponential term in 
Equation [15] essentially zero after several radiation length, making 



the photon spectrum independent of depth: 

7(e^,t) = 0.5677(e^,f = 0), (17) 

where the photon spectrum at the top of the atmosphere "f{e-,,t 
0) is as described in Equation 



14 



3 MUON PRODUCTION IN THE ATMOSPHERE 

High-energy 7-rays produce muons when they interact with the 
Earth's atmosphere. These muons will then traverse down to the 
bottom of the sea, producing Cherenkov light that can be detected 
by the detector array. This idea of detecting 7-induced showers by 
detecting the produced muons has been around for a long time. 
However, early calculations performed in the 1960s seem to indi- 
cate that 7-induced showers are muon-poor, having only less than 
10% the muon content of proton-induced showers jStanev, Gaisser[ 
'& Halzen 1985). These calculations are contradicted when muons 
were firmly detected at underground detectors, coming from the 
direction of Cygnus X-3 (e.g. [Marshak et al.|[T985^ . Despite the 
low rates and weak signals, these detections raised the interest to 
build large-area detectors that can detect high-energy muons and 
thus operate as 7-ray observatory. |Stanev, Vankov & H alzen ( 1985|l 
then identify two channels in which muons can be produced in 7 
showers: photoproduction and direct muon-pair production. In pho- 
toproduction, muons are produced from the (semi)leptonic decay 
of pions or kaons produced by the interaction of high-energy pho- 
tons with the atomic nucleus of the atmosphere. This is the most 
important channel to produce muons in the GeV regime. In direct 
muon-pair production, muons are created directly via the channel 
7 + ^ — )■ ^ + /i^ + in which Z is a nucleus of the atmo- 
sphere. Whereas muon production through photoproduction dies 
away with increasing energy, the cross section for muon-pair pro- 
duction increases with energy and thus muon-pair production is the 
dominant muon producing channel in the TeV regime. 

In the following subsections we will describe the neces- 
sary formulation to calculate the muon flux generated in gamma- 
induced showers. For convenience, all units of length are converted 
into radiation lengths in the air Arad, which is taken to be 37.1 g 

„™-2 



3.1 Pion decay 

The interaction of high-energy photons with atomic nuclei in the 
atmosphere can produce pions through the reaction 7+ — >■ tt+X 
followed by leptonic decay of pions into a positive muon and a 
muon neutrino, or a negative muon and a muon antineutrino : 



(18) 



with a probabilty of close to 100% to occur. The formulation to 
calculate the muon spectrum from this chaimel has been calcu- 
lated using the linear cascade equation and assuming a power-law 
photon spectrum with spectral index 6 = 1 by Dr ees, Halzen &[ 
Hikasa ( 19891, and its generalisation to an arbitrary spectral index 



by , Halzen, Kappes & 6 Murchadha 1 2009| 



For the case of b 7^ 1, this paper will closely follow that of 
Halzen, Kappes & O Murchadha ( 2009 1, which begins by an ansatz 



that the differential pion spectrum in the atmosphere can be factor- 
ized as 

7r(e,t) =7(e,i = 0)7r2(e,i), (19) 
in which 7r2 (e, t) can be split in two regimes: the high energy 
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regime where pion interactions dominate over decay, and the low 
energy regime where pion interactions are neglected. The pion 
spectrum at high energy is 



HE/,-, 



exp(Aii) — exp(— t/A^) exp(A2t) — exp(— t/A^) 
(ao + Ai)(Ai + ^) 

Z-yTT + Al)((TO + A2) 



A-yA A2 — Al 

while the spectrum at low energy is 



(o-o + A2)(A2 + ^) 

(20) 



LE/ ,x 

T2 (e,*) 



(ao + Ai)(cro + A2) 



X / dt' 
/() 



A2 — Al 

oxp(Ait') exp(A2t') 

_ (70 + Al (70 + A2 

in which 5 — t/d,r, where is the decay length 

et cos 6 



(21) 



(22) 



here e,^ = 115 GeV is the pion decay energy constant. 

The integral in Equation[2T|can be expanded into series: 



dt' 



t'y exp(Ait') 



tj ao + K 0-0 + A, ^ (j - + j) 



In Equation[20land[27] 

Att = 173 g cm~^ = 4.66 radiation lengths 
is the effective pion interaction length in the atmosphere, 



is the ratio between cross sections a- 



A, 



.7A 



and (T-yjv, and 
446.14 radiation lengths 



(23) 
(24) 

(25) 
(26) 



= 7(<;,i = 0) min [7r™(t), 7r2'^(e, t)j 



(27) 



The muon flux at the surface of the Earth can then be obtained by 
using standard 2-body decay kinematics, assuming no muon decay 
and energy loss in the atmosphere: 

lo J. (l-r)e' dn{t) 

in which r = {m^/rriT^)^ and _B^^ = 1 is the number of muons 
produced for each decaying pion. The maximum depth tmax is de- 
termined using 

nax (■^) 7 — ^/i 



dN^ 

d€a 



(29) 



is the interaction length of photons in atmospheric nuclei. These 
values are assumed to vary little for different spectral indices and 
energy. 

Due to the unavailability of an analytical expression for both 
energy regime, taking a smooth transition from one regime to an- 
other is difficult. The pion spectrum at all energy regime is then 



where A^+g- = 9/7 is the electromagnetic cascade length and 
(a:)-y_).^ = 0.25 is the fraction of 7-ray energy that goes into the 
final muon for the case of pion decays. 

For the special case of 6 = 1, we calculate the muon spectrum 
using the formulation by |Drees, Halzen & Hikasa] ( |1989^ : 

de^ XjA 1 + [L~f/H~, je^e^r cos 9 




Figure 3. Feynman diagram for lepton-pair production in the presence of a 
nucleu.s A' 



where 



2(1 -r) A^ 



1 



3(1 -r) 



1 + ln 



(31) 



The constant terms (A^, z^t,, X^a) in the Equations above are the 
same as in Equations |24j|26] 



3.2 Direct muon-pair production 

The Feynman diagram for direct lepton-pair production 7 + — > 
N + + is pictured in Figurels] This reaction occurs when an 
impacting photon interacts with a photon within the electric field of 
a nucleus, producing a pair of leptons. The second photon is nec- 
essary to maintain the conservation of 4-momentum, transferring 
the required momentum from the nucleus. Lepton-pair production 
is related to bremsstrahlung by a substitution rule and the calcula- 
tion of the cross section can be done if we know how to calculate 
bremsstrahlung by electrons ( |Tsai|1974l l. For the interaction of a 
photon with nuclear electrons to produce muon-pair, the photon 
energy threshold must then be 



2mu 



(m^ + me) ~ 43.9 GeV, 



(32) 



where rrie is the electron mass and is the muon mass. 

To calculate an approximate formula of muon-pair production, 
what is usually done is taking the Bethe-Heitler result for electron- 
pair production i jBethe & Heitler|1934^ and substitute the electron 
mass with that of muon. This generalization would not be correct, 
however, because the atomic form factor involved in the calculation 
must be integrated over the transferred momentum in which the up- 
per limit is approximately the m ass of the lepton involved ^Halzen^ 
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We will now discuss the necessary calculations to obtain the 
accurate formula for the cross section of muon-pair production. 

The impacting photon energy will be fully shared by the re- 
sulting muon-pair according to 

e7 = 4+<^M' (33) 
or in terms of fraction of photon energy: 

x+ = —, X- — —, x+ + X- — 1. (34) 

To take into account the atomic and nuclear form factors, we 
need the differential cross section equation as a function of the 
momentum transfer Since this work concerns very high-energy 
photons, we can use the ultrarelativistic approximation written as 



On the detection of TeV ^ -rays from GRB 7 



( |Bethe&Heitler|1934t 

da 

dX4 



ro- 



a; I + x_ 



+ -a;+a;_<l>2(5) 



(35) 

where alpha is the fine-structure constant, Z is the charge of the 
nucleus — for the Earth's atmosphere Z — 7.37 (Rossi' 1 932l l, ro is 
the classical electron radius, and S is the screening parameter equal 
to the necessary minimum momentum transfer from the nucleus: 

6 ~ Qmin = • (36) 

The functions $12 are integrals of form factors over trans- 
ferred momentum q. Whereas electron-pair production involves 
only the atomic form factors, in the case of muon-pair production 
it is also necessary to consider the nuclear form factors since the 
momentum involved is much larger than the inverse square of the 
atomic radius l |Tsai| 1974^ 1. The functions $1,2 would then be 

$i,2(<5)= r°"^[F4g)-K(Q)]^ 1/^1,2(9, 5), (37) 
Js 1 

where F„ and Fa are respectively the nuclear and atomic form fac- 
tors and '01,2 are the wave functions of the nucleus. 

Equation [37] has been solved with several assumptions. We 
take the solution of jKelner, Kokoulin & Petrukhin| ( [T995[ l in which 
a single function $(5) = $1 = $2 is used for the case of com- 
plete screening. By taking the effects of complete screening into 
account we consider the fact that atoms are essentially neutral at 
large distance. This is because the electric charge of the nucleus 
get "screened" by the atomic electrons, i.e. their field are canceled 
by opposite electric charge of the atomic electrons, reducing the ef- 
fective charge according to distance and thus limiting the maximum 
distance at which photons can still interact. 

The contribution from inelastic form factors is also con- 
sidered. This must also be taken into account since muon 
bremsstrahlung occurs on electrons bound in the atom and not on 
free electrons ( Keln er, Kokoulin & Petrukhin|1995l l. 

Having considered both elastic and inelastic form factors, 
Equation|35]then becomes 

2 



da 

dx 



ro- 



1 



-x{l 



$cl(5) + ^$in{5) 



The elastic contribution $ei(5) is in the form of 



$ci(<5) = In 



$0 



1 + (Dr, 



,1/2 



1 + BZ-i/3ei/25/me 



Z 
(38) 



(39) 



where 



$0 



BZ 



-1/3 



1/2 



B = 202.4 
B = 183 



D„ = 1.49 

i3„ = 1.54A°-2'' 



1.6187. 



for Hydrogen, and 
otherwise. 



(40) 

Here A is the atomic number of the nuclei involved. For our case 
of the Earth's atmosphere, A = 14.78 ( |Rossi|1952) . 
The inelastic contribution <&in((5) is 



In 



In 



B'Z-^/^e^/'^S/nie 



(41) 

where B' — 1429. We can see that the differential cross section is 
symmetric in x+ and a;_ , thus we can write 

J. ^ _ ^ _ ^2 




Figure 4. Differential cross section of muon-pair production (Equation[38j 
in the Earth's atmosphere for various values of impacting photon energy 
as a function, as a function of x = e^^/ej which is the ratio between the 
resulting muon energy and the photon energy. The atomic and mass number 
of the atmosphere is taken to be (A, Z) = (14.78, 7.37). 



where a; substitutes either x-i- ora;_ and the other becomes (1 — x). 

In Figure [4| Equation |38| for various values of photon energy 
is shown. We can see that due to the "screening" effect the cross 
section does not increase indefinitely but saturates as e-, increases. 
I integrate the differential cross section over x to obtain the total 
cross section as a function of photon energy and the result is shown 
in Figure [5] In the figure it is shown that saturation of the cross 
section occurs when the impacting photon energy e-, ~ 10 TeV. 

Using the cascade equation, we can calculate the muon-pair 
flux at sea level: 



2A,.ad-;^7o I 



dxx 



da 

dx 



dt-i2{t,b), 
(42) 



where A'^^ is the Avogadro number. 



3.3 Other channels of muon production 

A 7-shower can also produce kaons and the hadronic decay of 
kaons can produce a positive muon and a muon neutrino or a nega- 
tive muon and a muon antineutrino: 



(43) 



This reaction, however, has only ~ 63.5 % chance of occuring 



Halzen, Kappes & O 



([Gaisser"1 990f . Furthermore, results from 

[Murchadhal l |2009) showed that the muon yield from kaon decays 
and other channels involving kaons can be neglected. 
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S [TeV] 

Figure 5. Total cross section of the process 7 — > fJ-'^ tJ-~ in the Earth's 
atmosphere as a function of impacting photon energy . Due to screening 
effect which limits the maximum distance in which high-energy photons 
can still interact with the nucleus, the cross section saturates for impacting 
photon energy £7^10 TeV. 



Positrons produced in 7-showers can also produce pairs of 
muon by interaction witii an atomic electron through reaction 
e'^e^ — >■ + . However, cross section for this reaction is very 
small and peaked at --^ 61 GeV and f alls rapidly with energy and 
is essentially zero for > 700 GeV jflalzen, Kappes & 6 Mu~ 
|chadha|2009| l. Thus this production channel can also be neglected 
altogether. 



3.4 Cosmic ray-induced muon background 

In order to calculate the detection significance of photon-induced 
muons, we need to know the amount of the background in our 
observation. In our case of photon-induced muons detection, the 
background consists of cosmic-ray induced muons. These muons 
are produced mainly through leptonic decay of pions, which is es- 
sentially the same channel discussed in Section [3T| Leptonic decay 
of Kaons is also another channel of muon production albeit it is less 
important. 

The energy spectrum of cosmic -ray induced muons, as a func- 
tion energy and zenith distance, has already been parametrized by 
|Gaisser, ( ,1990^ as 



0.054 



1 + 



1 + 



(44) 



This parametrization overestimates the actual measured muon flux 
for energies below 10 GeV because at that energy regime muon de- 
cay and muon energy loss become important factors (see Figure 6.1 
in |Gaisser|1990^ . However, this will not be our concern since this is 
far below the energy regime we are interested in, and Equation |44| 
fits perfectly well for high-energy regime. This equation estimates 
the muon flux at sea level, thus if we want to estimate the muon 
background at detector we have to apply the appropriate muon en- 
ergy loss formula for seawater. We will discuss this later in Section 

H 



4 PASSAGE OF MUONS THROUGH SEAWATER 

Upon traversing a medium, energetic muons lose their energy 
through ionization and radiative processes. This energy loss can 
be treated by taking the standard formula to calculate the average 
energy loss jBarrett et al.|1952) 

- J = a{e) + h{e)e, (45) 

in which a(e) is the ionization contribution of the energy loss, while 
fe(e) = &p(e)+6i,(e)+6„(e) is the radiative contribution consisting 
of e^e~ pair production {bp), bremsstrahlung (6b), and photonu- 
clear interaction 6„ . 

Here I take the approach of [Klimushin, Bugaev & Sokalski| 



( 2001^ by splitting a(e) into two separate processes, a{t) 
ac(e) + ae(e), where Qc is the classical ionization process suffi- 
ciently described by the "Bethe" equation l |Nakamura & Particle| 
[Data Group 2010) and is the e diagrams for bremsstrahlung 
treated as part of an ionization process. Oc can thus approximated 
by 



= + flci In 



rua 



(46) 



1 + 



in which VKmax is the maximum transferable energy to the elec- 
tron and m^,e are respectively the masses of muon and elec- 
tron. The coefficients, in units of (10~® TeV cm^ g~^)> 
(accflcj = (2.106,0.0950) for e ^ 45 GeV and (0^0, acj = 
(2.163, 0.0853) for e > 45 GeV. For Be, a polynomial approxi- 
mation is used: 

aAe) = 3.54 + 3.7851ne + 1.151n^e 

39 21 (47) 

+ 0.06151n e 10"^ TeV cm 

where e is in units of GeV. 

The terms of b are parametrized in a polynomial function in 
the form 



bi (e) = bij W e, where i = p,b,n. 



(48) 



Here the energy input e is also in units of GeV. The values of co- 
efficients for bij is already calculated by [Klimushin, Bugaev &| 
|Sokalski| ( |2001| l and is tabulated in their Table II. These formu- 
lations of energy loss are expected to still valid for edotcctor ~ 
30 GeV — 5 TeV and slant depth (3 — 12) km with errors up to 
±(6 - 8)% (|KlimushrnTBugaev&~Sokalski|2001|l. 



Taking into account these contributions, the total muon energy 
loss in seawater as a function of energy is shown in Figure |6] In 
this figure we can see that at high energies radiative processes are 
more important than ionization. The critical energy at which the 
energy loss from ionization and radiative processes are equal can be 
calculated by solving e^c ~ cL{e^c) /6(eMc). In the case of seawater 
this is e^c ~ 590 GeV. Below this critical energy the dominant 
process is ionization while above this limit the radiative processes 
starts to dominate. 

If we integrate Equation|45]we can obtain the integral equation 

de 



a(e) + fe(e)e 



+ R = Q, 



(49) 



in which esurfaco is the energy at the surface of the sea and edotcctor 
is the energy at detector level, located at slant depth R = d/ cos 6 
where d is the vertical distance of the detector and 9 is the zenith 
distance from which the source came. Solving this equation, we 
can obtain esurfaco if edctcotor is the input and vice versa. I solve 
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Figure 6. The muon energy loss in seawater as a function of energy, calcu- 
lated from Equations |46] to |48] The total energy loss (solid line) is decom- 
posed into contributions from different processes, indicated in the legend. 
This Figure is made using the values of [Klimushin, Bugaev & Sokalski| 
(200T) . 
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Figure 7. The muon energy loss by passing a layer of sea water with vertical 
depth d = 2475 m is pictured here in the form of muon energy at the 
surface of the sea Egurfaco as a function of muon energy at the detector 
level edetector- We plot the energy loss for different zenith distance 6, thus 
the path length is R = d/ cos 9. 



Equation |49] to obtain esurfaco as a function of edetector- The result 
for ANTARES depth of d = 2475 m below sea level is shown in 
Figure|7]for several slant depths. 

The relation between tsurface as a function of edetector is par- 
ticularly useful to obtain the muon flux at detector level: 



5 MUON FLUX FROM SINGLE GRB 

Once we know how to produce gamma ray-induced muons in the 
atmosphere and how they lose their energy in seawater, we are now 
in the position to calculate the muon yield both on the surface of 
the sea and at detector level. I first calculate muons produced from 
a fictive, unattenuated test source with fluence f-, = 10"^ TeV"^ 
km~^ s~^ at 1 TeV. The source is a point source with negligi- 
ble diameter, assumed to be located at zenith distance 6 — 30°. 
The muon flux is calculated for three alternatives of spectral in- 
dices b = (0.6, 1, 1.6) and cutoff energy at emax ~ 300 TeV. For 
the background estimation, the opening angle of the search cone is 
taken to be Scone = 1°. The results are shown in Figure [I] and com- 
pared to a background of cosmic ray-induced muon flux calculated 
for the same zenith distance. 

These results are larg e ly con sistent with the results of |Halzen^ 
Kappes & O Murchadha 1 2009 1. We can see that the dominant 



dN 
dedet 



(edet, -R) = 



dN desui 



de 



dede 



(50) 



channel of muon-production at low energies is by pion decay. How- 
ever the number of muons that can be created from this channel 
goes down with photon energy. At high energies, because the cross- 
section of the muon-pair production goes up with photon energy 
before reaching saturation point at e-y > 10 TeV, the dominant 
muon production mechanism is direct-pair production. 

Confident with consistency of the calculation, I proceed by 
calculating the muon flux for single GRB events located at differ- 
ent redshifts. Using Equation [14] I calculate the photon flux ar- 
riving at the top of the atmosphere from GRBs with spectral in- 
dices & = (0.5, 1, 1.25, 1.5), redshifts 2 = (0.05,0.1,0.2,0.5), 
and zenith distances cosO = (0.5, 1). A typical GRB power spec- 
trum measured by BATSE is 6 ~ 1.25 jPreece et al.|2000 |, however 
measurement inconsistencies has been reported and thus the shape 
of the spectral index at high energy is still debatable and might not 
be in the form of a simple power law (see e.g. [Kaneko et al.|2008] 
and [Gonzalez et al.|2003| l. Until this debate is clarified by Fermi, 
it is reasonable also to asume a milder spectrum with index 6 ~ 1. 
The other spectral indices, the shallower b — 0.5 and the steeper 
b — 1.5, while not entirely impossible nevertheless have a small 
possibility of occuring and is thus also considered to study their 
possibility of observing the muon signal. 

Throughout the calculation, the values Ai, = 10 s, ebk» ~ 
(b - a)epk./(l - a) = (b - 1)400 keV, and LZ^, = 8.9 x 
10^^ erg is used. The values taken for At», epk*, and are all 
the mean values determined from Swift results jButler et al.|2007] 
[Butler, Bloom & Poznanski|2010 '). After calculating the number of 
photons at the top of the atmosphere, the muon flux at the surface 
of the sea is then determined by means of Equation |28| or [30] — 
depending on the spectral index considered — and Equation[42| The 
muon flux at the surface is then transformed to the muon flux at 
detector level by way of Equation|50| and the corresponding energy 
at detector level is calculated by solving Equation [49 [ 

The results of this series of calculations are shown in Figurel 
using the attenuation model by [Finke, Razzaque & Dermer[ ( [201()| 
One panel in each of these Figures plot the muon flux of GRBs 
for one spectral index. For each spectral index, the muon flux from 
GRBs at different redshifts is also shown and indicated with the 
color scheme shown in the legend. For each redshift, an area is 
drawn to show their dependence on zenith distance. The the bor- 
ders of the area drawn for each redshifts are the the muon flux at 
zenith distance 9 = (solid lines) and at 9 = 60° (dashed lines). 
Anything in between those two lines are then the amount of signals 
from any zenith distance between the borders. A background flux 
consisting of cosmic ray induced-muons calculated from Equation 




Figure 8. The vf^ spectrum of a fictive, unattenuated test source with fluence = 10^^ TeV^^ km^^ s^^ at 1 TeV. The spectrum is decomposed into 
its major contributing components: Pion decay and direct pair production. I calculate the spectrum for three different photon spectral index s = (0.6, 1, 1.6) 
and compare it with a background spectrum of cosmic ray-induced muons calculated using Equatio n | 44 | assuming that this fictive source is located at zenith 



distance d = 30°. The search cone has an opening angle of 1°. The result is largely consistent with Halzen, Kappes & O Murchadha 



2009 



|44|is also shown for the same limit of zenith distances. In the Fig- 
ure it is indicated by the black area. The search cone or the opening 
angle is taken to be 1° . I performed the same calculations for other 
attenuation models, but upon inspection of the numbers, results in- 
dicate that the magnitude of attenuation does not differ much for 
nearby universe, i.e. 2 < 0.2. Hence here only results calculated us- 
ing the calculation by lFinke, Razzaque & Dermer| ( |2010t is shown. 

For b < 1, we have a situation in which at low energy the 
muon yield at slanted angle is higher than the yield at vertical an- 
gle for the source located at the same distance. This is because ac- 
cording to the Cascade Equation discussed in Section [23] a pho- 
ton shower with a shallow spectrum — i.e. having spectral index 
b < 1 — will instead grows instead of dissipating in the atmosphere. 
Because at slanted depth the photon passed through thicker layers 
of atmosphere, more photons are created in the cascade and thus 
have more chance to produce muons. At higher energies that situa- 
tion is however no longer the case because the effectiveness of the 
pion decay channel in creating muons goes down with energy and 
the cross section of the muon-pair production channel reach satu- 
ration point. At slanted angle this limited number of high-energy 
muons then still have to pass a thicker layer of seawater than those 
that has to be traversed by muons coming from vertical direction. 

The results shown in Figure [9] indicate that the number of 
muons reaching detector level depends heavily on the GRB's dis- 
tance from us and its power spectrum. The redshift becomes the 
most important factor because it determines the number of photons 
that survives all the way from the GRB to the top of the atmosphere, 
the power spectrum comes second as it determines the number of 
photons originally produced in the GRB. 

The muon spectrum is then integrated to obtain a muon event 
rate with energies higher than c^t, detector- 

f°° dN 
N^{> e^.dctector) = / de^---^ (51) 

■/ej, , detector '^^M 

The result of this integration is shown in Figur e |10| u sing the at- 
tenuation model by |Finke, Razzaque & Dermer| ( |2010) . This result 
can give us an idea of how many muon events per unit area per 
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Figure 9. The muon flux at detector depth (2475 m) for single GRBs emit- 
ted from different redshifts as indicated by the color coding on the legend. 
The color black is the background flux from cosmic ray-induced muons, 
calculated using Equation |44| assuming a search cone with an opening an- 
gle of 1°. For each color, the muon flux drawn by the dashed-line is the 
flux from zenith distance 8 = 60° while those drawn by the solid line is 
the flux straight from the zenith (i.e. 8 = 0). The filled-area then defines 
all the possible flux from all zenith distance between 6 = and 8 = 60°. 
Attenuation is determined by using a model by |Finke, Razzaque & Dermer| 
j20l0j. 



unit time can we expect from any GRB event with the given power 
spectrum, redshift, and zenith distance. 

To explore further the effect of distance to the muon event rate 
at detector level, in Figure[TT|the event rate of muons with energies 




10"^ 10"' 1 10 10"^ 10"' 1 10 10^ 



e,.,deteclor [TeV] 

Figure 10. The expected number of muons with energies higher than the 
given muon energy e^^dctcctor at detector depth (2475 m) for single GRBs 
emitted from different redshift. The color black is the background flux. The 
same as in Figure|9] dashed-lines indicate muon count from GRBs occuring 
at zenith distance d = 60°, while solid lines are muon count from GRBs 
occuring at the zenith (i.e. 9 = 0). In this Figure, the EBL attenuation is 
calculated using the model by |Finke, Razzaque & Dernier I _2010j . 
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Figure 12. The total number of muons with energies e^, detector > 
0.1 TeV for GRBs from different redshifts and different spectral index 
as indicated by the color code in the legend. The intrinsic burst duration 
At* is assumed to be 10 sec, thus making tgo = (1 + 2)At,. The total 
muon count is calculated by assuming different detector size, which is also 
assumed to be independent of energy. 




Figure 11. This figure plots the muon count with energies ep, detector > 
O.lTeV for GRB sources occuring at different redshifts. As in Figure [5] 
dashed lines are for GRBs at zenith distance d = 60° while solid lines are 
for GRBs at S = 0. EBL attenuation is calculated using the model by |Finke,] 
[Razzaque & Dermer|j20T0) . 



higher than 0.1 TeV per unit area per unit time, A^/j (e^, detector > 
0.1 TeV), is plotted as a function of redshift. Here the black lines 
are the background rate from cosmic ray-induced muons at zenith 
distances 6 = [0, 60°] and it is independent of redshift. 

Figure[TT|tells us the minimum redshift and maximum zenith 
distance to observe, for example, at least 1 muon event per kilo- 
meter square per second. For example, a GRB event with power 
spectrum 6 = 1.25 that occurs at the zenith must have a redshift of 
z < 0.07 if we want to observe at least 1 muon event per kilometer 
square per seconds. The number of muons produced from a photon 
spectrum with b = 0.5 and those from 6=1.5 exhibit a large dis- 
crepancies, ranging from A'^ ~ 1 to A^^ ~ 10* km~^ s~^. This 
is because a photon flux with a shallow spectrum can produce an 
electromagnetic shower that grows in the atmosphere, while flux 
with a steeper spectrum produce a shower that instead dissipate in 
the atmosphere. 

The number of detectable muons depends also on the size of 
the Cherenkov detector. ANTARES is projected to have an effec- 
10~^ km^ whil e IceCube is expected 
1 km^ (Halzen, Kappes & 6 



live muon area of A[,^^ 
to have an area the size of A^„ 

ett 



|Murchadha|2009t . In Pi gure[T2] I calculate the total number of de- 
tectable muons during the whole duration of the burst for four dif- 
ferent detector size. The downgoing muon effective area is assumed 
to be yl^g — (10~^, 10~^, 0.1, 1) km^, and is also assumed to be 
constant with respect to the muon energy. One quarter of Figure [T2| 
plots the total number of events for different detector sizes. 

We can see that although larger detector with A'^g = 1 km'^ 
detect more unwanted background muons, they also see farther 
GRBs, up to z ~ 0.3 for b = 0.5. A detector the size of ANTARES, 
however, can only detect at least 1 muon event from a GRBs at red- 
shift up to 2 ~ 0.2 for the same integral index. 

Since we know the number of signals and noises in our de- 
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Figure 13. The muon signal detection significance for GRB sources located 
at different redshifts and different spectral index as indicated by the color 
code in the legend. The detection significance is calculated using the |Li &| 
|MaH1983) formula, and using the quantities shown in Figure [12] 
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Figure 14. The combination of redshift z and spectral intex b that gives a 
detection significance of 3o" (blue lines) and 5o" (red lines), for GRB pho- 
tons that came from zenith distances of 0° (solid lines) or 60° (dashed 
lines). 



lector, we can now calculate the expected detection significance 
of each individual GRB as a function of redshift. The significance 
S is calculated according to the procedure outlined by ' Li & Ma| 
( |1983) . The total signal A'^on is the number of muon events within 
the Scone = 1° search cone and during the ton = tgo time interval, 
while the total number of background A'^ofi is the number of muon 
backgrounds within the same search area but some amount of time 
toff before the GRB took place. The statistical significance S (the 
number of standard deviation above background) is determined us- 
ing the likelihood ratio method: 



S = V2\ Non In 



1 + a 



iVo: 



+ NoB In 



(1 + a) 



iVoff 
iVon + iVoff 



1/2 



(52) 



where a is the ratio a = ton/^off- The time toff to measure 
the background rate is taken to be 2 hours, i.e. toff — 7200 s, 
thus making a very low. The result of this calculation is shown in 
Figure[T3] again with four different panels each assuming different 
detector sizes. 

The Figure predicts the detection significance of observing 
GRBs with a certain power spectrum, zenith distance, and red- 
shift. We can also use this result to determine the maximum red- 
shift where a GRB has to occur if we want to have at least 3a or 5a 
detection significance. As an example, for an ANTARES-sized de- 
tector to detect a GRB signal with 5a significance, a GRB event at 
zenith must be closer than z < 0.05 if its power spectrum is fe = 1. 



6 CONCLUSIONS 

The two most important factors in detecting the TeV components of 
a GRB are its redshift and its spectral index. Its redshift determines 



the number of photons that survived to the top of the atmosphere, 
while the hardness of its spectrum determines whether the electro- 
magnetic spectrum grows or dissipate in the atmosphere. The com- 
bination of these two values that give a detection significance of 
3a or 5a over the background is prese nted in Figure |14[ A typical 
GRB has a spectral index b = 1-L25 ( |Preece et al.|2000|[Kaneko| 
|et al.|2008| >. For an ANTARES-type telescope, a typical GRB must 
then be located at least at redshift z < 0.05, while a larger tele- 
scope with a muon collecting area of Af = 1 km^ can see further 
up to 2; < 0.1 for a GRB with the same characteristic. 

A recent analysis of Fermi GRB data by |Zhang et al.| |20TTJ 
suggests that the peak of the distribution in b has shifted instead to 
b ~ 1.6, a much steeper slope than what was suggested by previ- 
ous observations. Consequently, the maximum redshift that permits 
a 3a detection is lower: Redshift z < 0.005 for an ANTARES-type 
telescope and z < 0.01 for a km^ neutrino telescope. In the anal- 
ysis of |Zhang et al.| ( (2"01 l| l, the peak distribution of integral index 
a is a ~ —0.1, which is not significantly different with previous 
results. 

The limitation pertaining to distance proves to be a great hin- 
drance to the detection of TeV 7-ray from GRBs, as there are not 
many GRBs with known redshift that took place at so close dis- 
tance. Recent analysis of 425 Swift GRBs found out that the red- 
shift distribution of GRBs is peaked at z ~ 1 jButler, Bloom &| 
[Poznans ki 20101. Within this data set, there are 144 GRBs with 
known redshift and 3 of them have z ^ 0.15. This corresponds 
roughly to a fraction of P{z ^ 0.15) ~ 7 x 10~^. 

From these results we can conclude that the secondary role of 
neutrino-telescopes as a gamma-ray telescope can only be played- 
out restrictively to the nearest GRB sources. As nearby GRBs tend 
to belong to a different population (i.e. short GRB) than the ones 
further away, further considerations must also be taken in view of 
the different luminosity and burst duration of this population. 

The expected rate of muon signals calculated in this paper has 
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not yet include the detection efficiency of tfie detector. To under- 
stand tfiis detector effect a furtfier study must be done, wliicfi re- 
quire a Monte Carlo simulation of the interaction between the muon 
signals and the detector. This is a work reserved for future studies. 
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